Investigating Neutron Polarizabilities through Compton Scattering on 3 He 
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We examine manifestations of neutron electromagnetic polarizabilities in coherent Compton scat- 
tering from the Helium-3 nucleus. We calculate 7 3 He elastic scattering observables using chiral 
perturbation theory to next-to-leading order (C(e 2 Q)). We find that the unpolarized differen- 
tial cross section can be used to measure neutron electric and magnetic polarizabilities, while two 
double-polarization observables are sensitive to different linear combinations of the four neutron 
spin polarizabilities. 
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The theory that describes the internal dynamics of 
the neutron is quantum chromodynamics (QCD). The 
neutron has zero charge, but higher electromagnetic mo- 
ments encode the strong-interaction dynamics which gov- 
erns its structure. These quantities therefore provide 
tests of our understanding of QCD. For example, an early 
success of the SU (3) quark picture was its prediction of 
magnetic moments, jl, for the neutron and other strongly- 
interacting particles (hadrons). Magnetic moments are a 
first-order response to an applied magnetic field. In this 
paper we will be concerned with electromagnetic polar- 
izabilities, which quantify the second-order response of 
a particle to electromagnetic fields. The two most ba- 
sic polarizabilities are the electric and magnetic ones, a 
and /?, which measure the ability of an applied electric 
or magnetic field to produce an induced dipole moment. 
The Hamiltonian for a neutral particle in applied electric 
and magnetic fields, E and B, is then: 
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where we have worked up to second order in E and B, and 
have not yet considered terms which involve derivatives 
of these fields. For a spin-half particle consideration of 
first-order derivatives allows four new structures which 
are second order in E and B [lj. They are: 
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The coefficients 71-74 are the "spin polarizabilities". 
This paper will argue that for the neutron, the most ba- 
sic and stable neutral hadron, a, /3, and 71-74 can be 
extracted from Compton scattering on 3 He. 

Polarizabilities such as those in Eqs. ([T|) and @ can be 
accessed in Compton scattering because the Hamiltonian 
(fT| yields an amplitude for Compton scattering from a 
neutron target of the form: 
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Here t\— te are invariants constructed out of the photon 
momenta and polarization vectors (e and e'), and, in the 
case of t^-te, the neutron spin, e.g. t\ = e'-e and £3 = ia- 
(e' x e). The Ais are Compton structure functions. The 
to 2 terms of A\ and A2 involve a and (3 Q, while the lu 3 
terms of A^-A^ depend on 71-74 in various combinations. 

For the proton, an expression similar to Eq. ^ but 
supplemented by the Thomson term, — j^e' ■ e, applies. 
The larger cross sections that result from the addition of 
this term lend themselves to low-energy measurements 
from which and (3^ can be extracted. A consider- 
able number of 7P experiments over the past decade had 
this as their goal 0] . A combined analysis of their low- 
energy differential cross section (DCS) data yields Q : 
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(12.1±l.l(stat.))±8;|(th.) x l(T 4 fm 3 , (4) 
(3.4±l.l(stat.))lo;l(th.) x l(T 4 fm 3 . (5) 



No elastic Compton scattering measurement of the 

(v) 

7- s has yet been performed, but they affect double- 
polarization observables. Of these, A 2 and A x are de- 
fined by taking the beam helicity to be along 2; then A z 
(A x ) is the difference between the DCS when the target 
is spin-polarized along +z (+x) and along — z (—x). For 
lu < TOjr the 7j- P ' ) s affect A z and A x because of interfer- 



ence between A^ a{p) ^ a a(p) 



Aq ' and Af in the expressions for 



these observables An experiment which exploits this 

(r>) (r>) 

interference to probe 71-74 has been proposed for the 
High-Intensity Gamma-ray Source (HI7S) at TUNL 0. 

However, neither polarized nor unpolarized Comp- 
ton scattering experiments can be directly performed 
on the neutron, since it is not a stable target. A va- 
riety of techniques have been proposed to extract a'™' 
and f3^ n \ including neutron scattering from the Coulomb 
field of 208 Pb and Compton scattering on the deuteron — 
both elastic and quasi-free. The most accurate num- 
bers come from the last technique and yield (in units 
of 10- 4 fm 3 ) 0: 

a (n)_0(n) = (g.8±3.6(stat)±2.2(mod.)t 2 ; 1 (sys)). (6) 
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These numbers represent a fascinating interplay of 
long-distance (r ~ l/m f ) and short-distance (r ~ 1/A) 
dynamics. The dominant piece of is due to the 
"cloud" of virtual pions that surrounds the neutron. 
But there are also significant contributions from short- 
distance physics — especially in [3^ n \ This interplay can 
be systematically computed in baryon chiral perturba- 
tion theory (xPT), a low-energy effective theory that en- 
codes the low-energy symmetries of QCD and the pat- 
tern of their breaking (see Ref . [H for a review) . Observ- 
ables in are computed in an expansion in powers 

of Q = p '™ lr , where A is the excitation energy of the 
lightest state not explicitly included in the theory. At 
0{e 2 Q) there are no contributions to the 7n amplitude 
from a short-distance operator. The prediction for 
the 711 amplitude comes from nucleon-pole, pion-pole, 
and one-pion-loop diagrams, with the latter capturing 
the dominant piece of the "pion cloud". This 0(e 2 Q) 
calculation yields the entire dependence of A^—Aq on 
photon energy and scattering angle up to corrections of 
O(f). The 0(lu 2 ) and 0(lu 3 ) non-pole pieces of A ( { l) ~ 
4 n) then give @,[1]: 
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5e l 9 } = 12.2 x 10- 4 fm 3 
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= 2 7 < n) = 4 7 < n) = -4 7 j n) = 4.4 x 10- 4 fm 4 .(8) 

(The 7} s can also be written in terms of qa, fm an d 
m,m) The contributions of short-distance physics to 
Eq. ([7]) are suppressed by one power of Q, and to Eq. ([5]) 
are suppressed by two powers of Q. In addition, x?T pre- 
dicts that a^ p \ (3^ p \ and the 7| p ' ) 's are the same as the 
corresponding neutron quantities — at this order. These 
0{e 2 Q) predictions of %PT agree with the numbers in 
Eqs. (|4]) and ([6|) within the experimental error bars. 

We now examine how the predictions of Eqs. ([7]) and 
([5]) can be tested in elastic 7 3 He scattering. The scatter- 
ing amplitude is written as 



M = {9 f \d\9i), 



(9) 



with and f) being the anti-symmetrized 

3 He wavef unctions. The results quoted in this letter have 
been calculated using a wavefunction obtained from the 
Idaho-N 3 LO chiral potential 0] together with the NNLO 
chiral 3N force [T3| • For reviews of xPT applied to nu- 
clear forces see Ref. [llj. Note, however, that aspects of 
this power-counting are still under discussion [12j. 

The operator O in Eq. ^ is the irreducible ampli- 
tude for elastic scattering of real photons from the NNN 
system, calculated in xPT up to 0(e 2 Q). This is next- 
to-leading order (NLO), a lower order than was used to 
obtain |\&), and so our calculation is chirally consistent 
only to NLO. At NLO O has a one-body part 



with 1B (a) being the 7N amplitude where the external 



photon interacts with nucleon 'a'. O 



IB, 



(supplemented 



by what turn out to be very small corrections for the 
boost from the 7N cm. frame to the 7NNN cm. frame) 
follows from Eq. ^ and can be found in Refs. [l3j ]. 
Meanwhile the two-body part of O is: 
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0^(1, 2) + O zb {2, 3) + 0^(3, 1 
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and it represents a sum of two-body mechanisms where 
the external photons interact with the pair '(a, b)'. At 
0(e 2 Q) this operator encodes the physics of two photons 
coupling to a single pion exchange inside the 3 He nu- 
cleus. (We do not have to include any irreducible three- 
body Compton mechanisms in our calculation because 
they appear at the earliest at Q(e 2 Q 3 ).) We use the ex- 
pression for 2B given in Ref. [13( . This incorporates the 
few-nucleon physics that corresponds to the pion-cloud 
dynamics which yields Eqs. ([7]) and ||5J). As such it must 
be included on an equal footing with the polarizability 
effects that are our focus. The resulting 2B gives a sig- 
nificant contribution to the 7d cross section, and is an im- 
portant piece of the xPT calculations that provide a good 
description of the extant 7d DCS data 0, 0, 0, . We 
now simplify Eq. §9§ to: 



M 



3(*/l2(0 ls (l) + ls (2))+0 2B (l,2)|* 4 ), (12) 



using the Faddeev decomposition of The structure of 
the calculation is then similar for the one- and two-body 
parts. We calculate M. on a partial- wave Jacobi basis. 
Convergence of the results with respect to the angular- 
momentum expansion was confirmed. For details on the 
calculational procedure see Ref. fl6| . 
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FIG. 1: Comparison of different c.m.-frame DCS calculations 
at 60 MeV (left panel) and 120 MeV (right panel). 
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The amplitude (fT2|) is now used to calculate observ- 
ables. In Fig. Q] we plot our 0(e 2 Q) xPT DCS predic- 
tions for coherent 7 3 He scattering. The two panels are 
for uj = 60 and 120 MeV. Both show three different DCS 
calculations — 0(e 2 ), I A (Impulse Approximation) and 
0(e 2 Q). The 0(e 2 ) calculation includes only the proton 
Thomson term, since that is the 7N amplitude in %PT 
at that order. The I A calculation is done up to 0(e 2 Q) 
but does not have any two-body contribution. As ex- 
pected, we see that there is a sizeable difference between 
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the I A and the 0(e 2 Q) DCS: the two-body currents are 
important and cannot be neglected. Also, we see that the 
difference between 0(e 2 ) and 0(e 2 Q) is very small at 60 
MeV — showing that xPT may converge well there — and 
gradually increases with energy. This is partly because 
the fractional effect of aS n > and /?' n ) increases with u>. 
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FIG. 2: The c.m.-frame C(e 2 Q) DCS at 80 MeV with varying 
Aa'"' (left panel) and A/3 (n) (right panel). 

To quantify this, in Fig. H we plot the 0{e 2 Q) DCS at 
80 MeV obtained when we add shifts, Ao» and A^ n \ 
to the 0(e 2 Q) values of the neutron electric and mag- 
netic polarizabilities ([7]). We take Aa'"' in the range 
(-4... 4) x 10~ 4 fm 3 and Af3^ between (-2... 6) x 
10~ 4 fm 3 . This allows us to assess the impact that one 
set of higher-order mechanisms has on our 0(e 2 Q) pre- 
dictions. Two features of Fig. [2] are particularly notable. 
First, sensitivity to (3^ n ' vanishes at 9 = 90° because aS n > 
and (3^ enter in the combination + (3^ cos 9. 
Thus, and (3^ can be extracted independently from 
the same experiment. Second, the absolute size of the 
shift in the DCS due to Aa<"' and A/?( n ) is roughly the 
same for all energies. This suggests that measurements 
could be done at to w 80 MeV, where the count rate is 
higher, and the contribution of higher-order terms in the 
chiral expansion should be smaller. 

We have estimated the uncertainty due to short- 
distance physics in the three-nucleon system by using a 
variety of 3 He wave functions generated using various NN 
interactions with and without a corresponding 3N force. 
This produced changes of < 15% in the DCS at 120 MeV. 

Before examining double-polarization observables in 
7 3 He scattering we try to develop some intuition for the 
7 3 He amplitude. Since 3 He is a spin- 1 ; target the matrix 
element (|12l) can be decomposed in the same fashion as 
was the neutron's Compton matrix element in Eq. ([3]). 
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'He 



A\ B + Aj B , (13) 



where A\ B (A 2B ) comes from considering the matrix el- 
ement of the one-body (two-body) operators in Eq. (fT2|) . 
and the structures t^-t^ now involve the nuclear — not the 
neutron — spin. However, in 3 He the two proton spins 



are — to a good approximation — anti-aligned, so the nu- 
clear spin is largely carried by the unpaired neutron [17 1 . 
We find that the 0(e 2 Q) two-body currents A\ B and 
A 2B are numerically sizeable, but A 2B -Aq B are negligi- 
ble. Hence, to the extent that polarized 3 He is an effec- 
tive neutron, we expect A i He — A^ for i = 3-6. Using 
Eq. (fT3")) to translate this into predictions for A 2 and A x 
shows that the effects of j[ n ^ '-74™'' will be enhanced in 
these observables by interference with A 1 He . But A 1 He 
is — at least at co rj 80 MeV — dominated by the contri- 
bution of the two protons, and so we anticipate a more 
marked signal from the neutron spin polarizabilities than 
is predicted for the corresponding 7d observables [l8j |. 

We emphasize that these arguments are meant only 
as a guide to the physics of our exact 0(e 2 Q) calcula- 
tion. Our 3 He wave function is obtained by solving the 
Faddeev equations with NN and 3N potentials derived 
from xPT. All of the effects due to neutron depolariza- 
tion and the spin-dependent pieces of 2B are included 
in our calculation of the amplitude © . This yields the 
results for A z and A^ shown in Figs. [3] and |U There 
we have proceeded analogously to our computations of 
the 7 3 He DCS, this time varying the neutron spin po- 
larizabilities and seeing the effect on A z and A x . Fig. [3] 
indicates that A z is quite sensitive to j[ n \ 72™ \ and 

(n) 

74 . With the expected photon flux at an upgraded 
HI7S such effects can be measured |19j . If this can be 
done as a function of 9 we can extract the combination 
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274™') cos 9. Turning to A^, Fig. U] shows 



that varying 7} or 74 produces appreciable effects in 
Aa: — but in a different combination to the sensitivity in 
A 2 . Use of different 3 He wave functions alters these pre- 
dictions for A x and A z by < 7.5%. For a more detailed 
discussion see [l6|. Thus, A z and A^ are sensitive to 
two different linear combinations of 7^ , 73 , and 74 
and their measurement should provide an unambiguous 
extraction of 7}"^ , as well as constraints on 72™^ and 74 . 

These 7 3 He scattering calculations are the first calcu- 
lations for this reaction, and there is significant scope 
for improvement. Computing of the NNLO (0(e 2 Q 2 )) 
pieces of the 7NNN operator O would allow a more de- 
tailed assessment of the pattern of convergence. When 
this is done we anticipate three kinematic domains where 
convergence may be slow. First, since we use the heavy- 
baryon formulation of xPT, the pion-production thresh- 
old is at w = m^, rather than in the correct position 
for 7 3 He scattering which is ^4 MeV lower. An esti- 
mate of the impact of this discrepancy on observables 
suggests a ~ 5% difference in the DCS, and < 3% in 
A 2 and A^, at 100 MeV. Second, the power counting 
we used is not valid at energies < ^jf-. For instance, it 
does not reproduce the correct uj = (Thomson) limit 
for the nuclear target, since the terms in O that restore 
that limit are higher-order effects when u) ~ m^. (A re- 
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FIG. 3: A z at lo = 120 MeV with (left-to-right) -f[ n) , j ( 2 n) , and varied one at a time. For 0(e 2 Q) j\ n>, s see Eq. © 
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FIG. 4: A a: (cm. frame) at lo =120 MeV when -y[ n) (left) and 
74"' (right) are varied one at a time. Legend as in Fig. [3] 



cent computation for 70! scattering verifies that they are 
indeed small for lo > 80 MeV [H|].) We therefore expect 
that assessment of these two classes of corrections, while 
an important check on our results, will not significantly 
alter them. We believe that the most important correc- 
tion will come from the inclusion of A(1232) degree of 
freedom. 7d scattering calculations which included such 
effects found a sizeable impact on the DCS at backward 



14] 



angles for w»100 MeV 

Our results for 7 3 He scattering are obtained from xPT 
NN & NNN interactions and 7N & 7NN operators, and 
are accurate to NLO in the chiral expansion. These first 
results on this reaction suggest that and (3^ can be 
extracted from the 7 3 He DCS and compared to results 
from 7d experiments at MAXLab [2(J. Meanwhile, two 
different linear combinations of 7^ , 72 > an d 74 can 
be constrained by measurements of double-polarization 
observables at facilities such as HI7S. This would provide 
new information on neutron polarizabilities. 
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